September 29, 1997 
To appear in ApJ on May 1, 1998 



as 
as 



(N 
> 



Implications of Cosmic Microwave Background 
Anisotropies for Large Scale Variations in Hubble's Constant 

Yun Wang, David N. Spergel, and Edwin L. Turner 



D I Princeton University Observatory 



Peyton Hall, Princeton, NJ 08544 
email: ywang, dns, elt@ astro. princeton. edu 



O ■ Abstract 
00 

o 

^ \ Low amplitude (linear regime) cosmic density fluctuations lead to 

spatial variations in the locally measurable value of Hq (denoted as Hi) 

O 

^\ = {Hl — Hq)/Hq, which are of order 3 — 6% (95% confldence interval) in a 

' sphere of 200 h~^Mpc in diameter, and of order 1 — 2% in a sphere of 400 h'^Mpc 



in diameter, for three currently viable structure formation models (tilted CDM, 



■ ACDM, and open CDM) as normalized by the 4 year COBE DMR data. 

However, the true matter distribution power spectrum may differ from any 
of the currently viable models. For example, it may contain sharp features 
which have escaped detection so far. The measured CMB dipole velocity 
(the Galaxy's peculiar velocity with respect to the CMB rest frame) provides 
additional constraints on the probability distribution of 6h that supplement 
our limited knowledge of the power spectrum. For a matter power spectrum 
which consists of the smooth power spectrum of a viable cosmological model 
plus a delta-function bump, we flnd that given the CMB dipole velocity, the 



- 2 - 



95% CL upper limit of \5h\ increases approximately by a factor of two, but the 
probability distribution of 6h is non-Gaussian, with increased probability at 
small Sh compared to Gaussian. Abandoning model power spectra entirely, we 
find that the observed CMB dipole velocity alone provides a very robust hmit. 



y{Sff)R < 10.5h~^Mpc/R at 95% CL, in a sphere of radius R, for an arbitrary 

power spectrum. 

Thus, variations between currently available local measures of Hq and its 
true global value of a few to several percent are to be expected and differences 
as large as 10% are possible based on our current knowledge of the CMB 
anisotropics. 
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1. Introduction 



In the standard gravitational instability scenario ( [Peebles 1980 , Peebles 1993| ) for 



cosmic structure formation, linear growth of density fluctuations is produced by, and 
produces, spatial variations of the expansion rate. In this rather generic scenario, such 
variations are the inevitable implication of the extremely large scale structures which have 



been detected in the galaxy distribution, primarily via redshift surveys (e.g., [DaCosta et al. 
ILauer fc Postman 1994 [Lin et. al. 199^ , fl'adros fc Efstathiou 1996| ). 



The connection of such expansion rate variations, which are often thought of as 



peculiar velocity fields, to the large scale density distribution is of direct interest ([Strauss fc 
[Willick 1995[ ) but is also potentially important as a source of systematic error in attempts 
to measure Bubble's constant Hq, the overall mean expansion rate of the Universe. In 
particular, if the expansion rate is correctly measured in a local volume which is not 
sufficiently large compared to the biggest significant cosmic structures, the strong possibility 
of a difference between it and the true cosmic value of Hq must be considered. 

This potential problem has long been known in principle, was emphasized in the 
context of modern structure formation models by Turner, Ostriker, & Cen (1992) and 
has subsequently been considered by several authors ( [Wu et al. 1995[ , [Nakamura et al.[ 



T995i [Nakao et al. T995| , |Shi, Widrow, fc Dursi 1996| , [Wu et al. 199^ , |Shi fc Turner"T997[ ) . 



The issue is of greater interest than ever at present for two reasons: First, measurements of 
Hq by conventional distance ladder techniques have improved dramatically and now achieve 
credible precisions of order 10 to 20 percent. It follows that a systematic difference between 
local and global expansion rates of a comparable fractional size are important sources of 
error. Second, there is suggestive, though not yet compelling, evidence that direct physical 
methods for measuring Hq on extremely large scales (out to redshifts of order unity) 
( Pirkinshaw, Hughes, fc Arnaud 1991[ , [Jones et al. 1993[ , [Schechter et al. 1997[ , [Kundic el 
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al. 1997) give a somewhat smaller value than the best distance ladder determinations which 
apply out to redshifts of a tenth or less ([Graham et al. 1997| , [Tonry et al. 1997| , [Eastman^ 



pchmidt fc Kirshner 19961) 



In this paper we attempt to exploit our best and most direct empirical information 
about large scale (linear) cosmic density fluctuations, namely observed anisotropies in the 
cosmic microwave background (CMB), in order to predict and/or constrain the expected 
resulting variations in the expansion rate. Specifically, we study these intrinsic fluctuations, 
5h = {Hl — Hq)/Hq, where is the local measure of Hq, in terms of the matter distribution 
power spectrum P{k). (Note that although galaxies are biased tracers of the mass density 
field, on the large scales which are of interest to this paper — 100 Mpc or larger, linear 
biasing is a reasonable assumption for conventional structure formation scenarios.) 

In fact, we present here three separate calculations of this sort, proceeding from the 
one which gives the strongest result but which is most model dependent to the one which 
is weakest but most robust (model independent). First, we simply base our input P{k) 
on standard structure formation models in terms of its shape and use CMB observations 
only to set the overall normalization or amplitude. Here, our results agree with recent 
work by Shi & Turner (1997). Second, we consider an input P{k) with a shape and 
amplitude constrained by CMB fluctuations; for some scales {k values) P{k) is effectively 
directly observed, but we allow for the possibility of strong features on those scales not 
directly sampled by currently available data by imposing limits based on the CMB dipole 
(the Galaxy's peculiar velocity with respect to the CMB rest frame). Third, we impose 
no constraints at all on the input P{k) (which does not even appear explicitly in this 
calculation) other than that it be Gaussian and not violate the aforesaid CMB dipole 
constraint. 

Section 2 contains general expressions for 6h and related variables. In Section 3, we 
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compute Sh for matter power spectra given by three viable cosmological models which 
are normalized by the 4 year COBE DMR data and satisfy constraints from large scale 
structure data. In Section 4, we study the effects of unexpected features in P{k) which may 
boost 6h by adding a delta-function bump to a smooth matter power spectrum given by a 
viable cosmological model. In Section 5, we apply Bayesian statistics to derive robust upper 
limits on 6h and related variables, using the CMB dipole velocity of v=627 km/s ([Kogut e"t 



al. 1993| , [b'ixsen et al. 191MD ; these upper limits are independent of the actual form of the 
matter power spectrum. Section 6 contains discussions and a summary. 



2. General Expressions 

For an observer at ri who measures the Bubble's constant by summing over the radial 
recession velocity divided by distance of N objects located at rj, (j = 1, 2, N) ([Turner 
Qstriker, fc Cen 1992| ) 



. / N _ Hjrj) - gp _ 1 ^ V, ■ (r, - r,) 



To find 6h for a sphere of radius R centered around x, we write ( [Shi, Widrow, fc Dursi 
[I996[) 

5^(x)« = [d'y^- W{y - X), (2) 

J -no |y - 

where W^(y — x) is the top hat window function with radius R, 

(47ri?3/3)-S |y-x|<i? 

w^(y-x) = <i (3) 

0, |y-x|>i? 
In linear perturbation theory (Peebles 1993, 1980), the Fourier component of 



v(y) = (27r) 3 / d^k exp(-?k ■ y) i 



IS 



00.60 TT 

^, = \^5A (4) 

IK 
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for an open or flat universe. 5k is the Fourier component of the density fluctuation. Hence 



^0.60 



ikx 



where 



(27r) 

3 



(5) 



Cix) = — sin a; — / dy 
X"^ \ Jo 



sm y 

y . 



(6) 



Note that jC{x — > 0) = —1/3, and jC{x) < 0, as expected; overdensities {Sk > 0) lead to 
infall, which leads to the underestimate of the Bubble's constant {Sh < 0). 



5h is a Gaussian random variable, with mean and variance 

('^^)« = ^i, dkP{k)[{kR)£{kR)f. (7) 
The power spectrum P{k) — The variance of density fluctuations in a sphere of radius 
R is 



.P. 



R 



27r2i?2 



dkP{k)l{kR) 



3n{kR) 



n ^ 2 



kR 



(8) 



Clearly, the fluctuations in the measured expansion rate stem directly from the fluctuations 
in matter density. 

We can define Ql = SnGpL/{3Hl). To lowest order in 5p/p and Sh, 
Sn = {flL - no)/no = 6p/p - 26h. We find 



{S'n)R 



dkP{k) I (kR) 



^3i{kR) 
kR 



n ^ 2 



27r2i?2 7o 

Note that the fluctuations in il^ are contributed additively by 5p/ p and 5h {^{x) < 0). 

= SnGpL/iSHl) is related but not identical to the locally measured Q (by counting 
matter or dynamical methods). 

The variance of pecuhar velocity v and bulk flow velocity Vr (variance of pecuhar 
velocity in a sphere of radius R) also depend on P(k) , 



(9) 



27r2 



^ / dkP(k), 

Jo 



R 



27r2 



dkP{k) 



kR 



(10) 
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Note that {{5p/ pY)r, {S'^r, and {'v'^)r are all proportional to integrals over 

P{k) multiplied by a window function W{kR)\ Fig.l shows these window functions. Note 
that the window function of {Sfj)^ is somewhat similar to that of {{Sp/pY)r and {Sq)r, 
but very different from that of (v|.). Also, {{6p/ pY)r, and {6q)r all decrease much 

faster with R (the radius of the observed volume) than (v|j,) (note that we have factored 
l/i?^ out of the integrals for (5|^)_r, {{6p/ p)^) r, and {Sq)r), while (v^) is independent of the 
size of the observed volume. 

Even though 6h is caused by the presence of non-zero peculiar velocities [see Eq.(l)], 
it has only weak correlations with the measures of the peculiar velocity field. 



3. Theoretical Power Spectra with CMB Normalization 

Let us consider cosmological structure formation models which simultaneously satisfy 
constraints from the observed LCRS power spectrum ( [Lin et. al. 1996| ), the Bubble's 
constant range of 0.5 ^ h < 0.8, cluster abundance results, and the reasonable assumption 
that LCRS galaxies are approximately unbiased on large scales relative to the mass 
normalization provided by the 4 year COBE DMR data. Following Lin et. al. 1996 , we 
assume that in the linear regime 



2 n'^^'^ 1 n^^ 

Pgalaxy,redshiftspace{k) — b j 1 + — ^ h — — j Pmass,real space,linear 

{k). (11) 



Three viable models are: (1) TCDM: = I, n = 0.7, h = 0.5, = 0.05, b = 1.3; (2) 
ACDM: Qo = 0.5 = Q^, n = 1, h = 0.5, Qbh^ = 0.015, b = 0.9; (3) OCDM: Qq = 0.5, n = 1, 
h = 0.65, rihh^ = 0.015, b = 0.9. 



The power spectrum of a model is given by 

m ^ (_ 



"^^'^ ^of-ArV^g), (12) 
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where Pq = 27i^ {10^ Shor ■ 90f ■ 3000"" ^ The 4 year COBE DMR data give ( [Bunn fc White 



19971) 



1 94 f^-O- 785-0.05 InCo 



1 95 ^2-0-35-0. 19 Inf^o 



exp[-0.95 (n - 1) - 0.169 (n - 1)^], Q = 1; 
n < l,n = 1. 

T{q) is the CDM transfer function, given by ( Pardeen et. ah 1986|) 

ln(l + 2.34g) p. , , . .2 , , .... ^4l-V4 



(13) 



2^(9) = "-^^ ' "•"-'^^ 1 + 3.89 q + (16. 1 g)2 + (5.46 q)^ + (6.71 g)^ 
2.34 g L 

where g = k/{hT), T is the shape parameter. We use ( |Hu fc Sugiyama 1996| ) 

r = noha^/'^e^^^, 

where 62.7 = 2.12QK/2.7K, and 



(14) 



(15) 



1 "2 ! 



ai = (46.9^]o/^')°■''° 

: (12.0 nohY''^ 



a2 



l + (32.1fio/i')'°-''' 
l + (45.0f]o/^2)-°-^«2 



(16) 



For the three models considered, T = 0.467 (TCDM), 0.223 (ACDM), and 0.299 (OCDM) 
respectively. 

Figs. 2-5 show {{6p/p)^)r, (5^)/?, and {v^)r, as functions of R, for the three 

models (TCDM, ACDM, and OCDM). Figure 2 indicates, for example, that a perfect 
measurement of the local expansion rate in a sphere of diameter 20,000 km/s (200 /I'^Mpc) 
would provide a 95 percent confidence interval determination of the global value of Hq with 
a width of 3 to 6 percent, depending on the adopted structure formation model. 



4. Smooth Power Spectra with a Feature on Unobserved Scales 



The COBE observations of CMB fluctuations on scales larger than 7° ( Bennett et 
al. 1996| ) directly probe the potential fluctuations on scales larger than 300 Mpc. It 
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is possible to smoothly connect the COBE measurements to observations of large scale 
structure. In fact, there is a family of cold dark matter models (tilted CDM, open CDM, 
Lambda-dominated CDM, mixed dark matter) that are consistent with both measurements 
of large-scale structure on scales of ~ 1 — 60 Mpc and CMB observations on large scales 
( IRatra et al.T997i Punn fc White 1997| ). 



There are, however, a number of observations that suggest that the amplitude of 
fluctuations on the 60 — 300 h^^ Mpc scale exceeds that predicted by CDM models or 
by smooth interpolations between the scales probed by galaxy surveys and the COBE 
observations: Lauer & Postman (1994) measure much larger bulk flows than predicted 
by standard cosmological models ( [Strauss et al. 19"^ ; deep pencil beam surveys detect 



excess power on scale of ~ lOO/Zi Mpc ( [Broadhurst et al. 1990| , [Cohen et al. 199"6D ; 



two-dimensional measures of the power spectrum in the Las Campanes redshift survey also 
find excess power on the lOO/Zi Mpc scale ( [Landy et al. 1996[) ; and K-band galaxy counts 
can be interpreted as indicating a very large scale local underdensity region ([Phillips fc 



Turner 1998|) . If the ratio of the baryon density to the closure density, Qbh^, is close to 



the low value suggested by the D/H measurements of Songaila, Wampler & Cowie (1997), 
then recent CMB measurements on the 1-2 degree scale ( [INetterheld et al. 1997 ) would also 



imply excess power near the 100/ h Mpc scale. Shi et al. (1996) have used the degree scale 
CMB measurement to constrain a possible feature in P{k) beyond about 100 Mpc scale; 
however, even if Q^h'^ is large, it is difficult to use CMB observations to place an upper 
limit on density fluctuations on this scale as CMB fluctuations on scales smaller than that 
probed by COBE may have been suppressed by reionization of the intergalactic medium at 
z > 20. 

These observations suggest that P{k) is not smooth, but may have some feature on 
the 60 — 300 Mpc scale. There are several interesting physical processes acting on these 
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scales because of the horizon size at the matter-radiation transition and at baryon-photon 
decoupling, so it is quite possible that there is new physics acting on these scales. Since 
we are trying to constrain the uncertainties in the Bubble's constant due to large scale 
structure, we will study the possibility of excess power by adding a delta-function bump to 
a smooth matter power spectrum Po{k) given by a viable cosmological model. Let us write 

P{k) = Po{k) + A5ik-ko), (17) 

We take Po(^) given by the tilted CDM model (TCDM) as a convenient smooth functional 
form, with n = 0.7, fio=l, = 0.05, and h = 0.5. Since dk Po{k) = 828.808 /^-2Mpc^ 
and A = %^^H^/{2tt^) [J;^ dk Po{k) + A] is less than (1831 km/s)^ at 95% CL (see the 
following section for more details), we find that A/ (271^) < 293.3 /I'^Mpc^ at 95% CL. 
Calculations show that the delta function spike has the most significant effect on the bulk 



flow velocity y{v'^)R- 

Now we compute the probability distribution of Sh measured in a sphere of radius R, 
given the CMB dipole velocity of f=627 km/s. Using the Bayesian theorem, we flnd 

P{Sh\v)r = ^ PiSH\A)RPiA\v) oc ^ PiSH\A)RPiv\A)PiA), (18) 
A A 



P{6h\A)r = ^ exp f-||^) . (19) 



where 



{^h)r is a complicated function of A. For the power spectrum of Eq.(|T^, we have 

(v^) = A = ai + a2A, 
{51) R = f],{R)+(32{R,ko)A. (20) 

We have defined 

01.20 0-2 .oo O1-20/T2 

27r2 Jo ^ ' 2n^ ' 
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/32 = ^[(^oit:)>C(A;oi?)]^. (21) 



The parameter A characterizes the departure of the power spectrum from the 

underlying smooth power spectrum; we expect the probabihty of A to decrease with 
increasing A. Let us use the prior probabihty 

1/A, A > Ae; 



P{A) = 



(22) 

0, A< A,. 



Ac is a cut-off motivated by physical considerations. We consider two choices for the cut-off 
Ac. 



(1) Ae = «i/a2 = /o°°dA;Po(A;); 

(2) ^ = 0.1 ai/a2 = 0.1 /o~ dA; P^ik). 

The first choice of Ac indicates the reasonable upper limit of 74c, while the second choice 
indicates a typical value of Ac of interest. 



The probability distribution of 5h given the value of v is 

= jfTm r (x+i)3/4+xov^ {- 



3f ^ X 6jjX 



2ai{x + l) 2pi{x + xi) 



(23) 



where 



..^"-f, (24) 

a2 Ac Q!2Pl 



and 

r^<^ X^l"^ { ?)V^ X 







The variance of bii given v — ^TI km/s is 

^^^'"^^ = ^ r xV2(J+' 1)3/2 {-2j(x + l)} • ^26) 
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Figs.6-8 show P{Sh\v)b^ for ko = 0.01 /iMpc"^ and R = 40, 100, 500 /I'^Mpc 
respectively. The sohd and dashed hnes are the distributions given by Eg. (]23|) , with Xc = 10 
{Ac = 0.101/0:2) and Xc = I {Ac = 01/02) respectively; the dotted lines are Gaussian 
distributions with the same variance [given by Eq. (p6D ] for Xc = 10. Note that P{Sh\v)r 
becomes increasingly non-Gaussian for increasing R. Table 1 lists the 95% CL upper limits 
on \Sh\, as well as 2a = 2^JJsj^\v)ji, the 95% CL upper limit on Sh if its distribution 
were Gaussian. 





R 


40Mpc/i-^ 


lOOMpc/i-i 


500Mpc/i-i 


Xc = 10 


2a 


0.115 
0.120 


5.50 X 10-2 
6.13 X 10-2 


9.44 X 10-3 
1.05 X 10-2 


Xc = l 


S'h 
2a 


0.135 
0.143 


8.67 X 10-2 
9.56 X 10-2 


1.59 X 10-2 
1.75 X 10-2 



It is not surprising that our results depend on the prior probability P{A) via the cut-off 
Ac, since different choices of Ac represent different input of physical information. It is worth 
noting that the 95% CL upper limit on \6h\ is close to 10% for R > lOOMpc/i-^ with 
the largest reasonable Ac (when the contribution to the mean square peculiar velocity by 
the bump in P{k) is greater or equal to the contribution by the underlying smooth power 
spectrum) . 

Also note that in our Bayesian statistics, the measured CMB dipole velocity of 
V = 627km/s has the effect of distorting the distribution of 6h away from Gaussian by 
increasing the probability of smaller 6h- 



5. A Robust Upper Limit on the Variance of 5h 



The CMB dipole velocity of f = |v|=627 km/s is the Galaxy's peculiar velocity 
with respect to the CMB rest frame. Note that v has three components which are three 
independent Gaussian random variables. In the previous section, it provided us with a 
constraint on possible unobserved features in P{k), but we can in fact use it in a more 
fundamental way to provide a single integral constraint on P{k) which is related in a simple 
way to possible expansion rate variations, as described below. We can infer the distribution 
of A = (v^) by using the Bayesian theorem. 




(27) 



P(A) is the prior probability of A. Since A/3 is the variance of the Gaussian variable 
{x component of v), it is most reasonable to choose -P(A) = 1/A. For diagnostic on the 



dependence on the prior, let us write 





A^^^ < aJ^^ at 95% CL, with aJ''^ given by 




(29) 



Using 




(30) 



we find that aJ^^ = 1831 km/s for A^^ < 200 km/s, and Aq^ = 2183 km/s for 

-1/2 ~ 1/2 

A^ = V = 627 km/s. Since Aq changes by less than 20% for the significant cut-off of 
A^^ = V, we can take A^^ = 0. 




/-oo dk P{k) 



27r2 Jo hU^c-^ /i-3Mpc 



3 < 



1831 km/s^ 
100 km/s j 



335.3 



(31) 
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at 95% CL. 

Using Eqs.(0), (g), and (g), we obtain at 95% CL 

{6UR)) < 



R J ' 

^SUR)) < .-o(l^^y. (32) 

We have used [xC{x)]^ < 0.3282, x'^{3ji{x)/xy < 1.7123, and x^[3ji{x)/x - 2/:(x)]2 < 
5.7579. 

This calculation provides us with the most robust results, both in the sense of being 
model independent and of relying only on the largest and most observationally secure CMB 
anisotropy (the dipole). They are thus also the weakest quantitatively. For example, 10% 
variations in the expansion rate are allowed on scales of 20,000 km/s (200 /i~^Mpc) in 
diameter at 95 percent confidence. 



6. Summary 

Cosmologists attempt to derive properties of the large-scale structure of space-time 
from local observations. These extrapolations rely on the assumption that we are probing 
a fair sample of the universe, so that physical quantities measured locally such as the 
Bubble's constant, or the mass-to-light ratio are representative of the universe as a whole. 

Large-scale structure generates deviations from the Hubble flow; thus, it is important 
that measurements of the Bubble's constant probe a large enough volume so that the 
effects of peculiar motions are small (|'l'urner, Ustriker, fc Cen 1992|). In this paper, we have 



estimated the expected variance in the Bubble's constant due to large scale structure. We 
began by considering currently fashionable models, which are consistent with galaxy surveys 
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on small scales and CMB observations on large scales. In these models, measurements of 
the Hubble's constant that are based on galaxies in a sphere of diameter 200 h^^Mpc are 
likely to be within 3 — 6% (95% confidence interval) of the global value, and of order 1 — 2% 
in a sphere of 400 h~^Mpc diameter. These limits assume that the current set of large-scale 
structure models are good approximation to the primordial power spectrum. 

The CMB dipole velocity (the velocity of the Galaxy with respect to the CMB rest 
frame) places a strong constraint on the amplitude of large scale structure. If there were 
enormous local void and density fluctuations, as suggested by several authors ([Harrison 



1993| , |Wu et al. 1995| , |Wu et al. 1996| ), then we would expect that the Galaxy would be 



moving with respect to the microwave background rest frame at a velocity much larger 
than the meaured value of 627 ± 22km/s ( P:<.ogut et al. 199^ , [I''ixsen et al. 199^ . Thus, 



the measured CMB dipole velocity can be used to derive strong constraints on density 
fluctuations on scales larger than those probed by redshift surveys. We have used these 
constraints to calculate the variations in the Hubble's constant, 6h, for power spectra which 
have a sharp bump on unobserved scales; we find that the 95% CL upper limit of \6h\ 
increases approximately by a factor of two. 

Finally, we have used the constraints derived from the CMB dipole velocity to place a 
robust limit on variations in the Hubble's constant. With the CMB dipole measurement 
alone, we are able to constrain variations in the Hubble's constant to be less than 10% 
on scales of 20,000 km/s (200/i~^Mpc) in diameter at 95 percent confidence. Thus, 
measurements of Hq that probe out to this scale are likely to be accurately probing the 
global expansion rate of the universe. 
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Fig. 1. — Window functions of {Sjj)ji (solid line), {{Sp/p)'^)ji (dotted line), {Sq)r (short- 
dashed line), and (v^)^ (long-dashed line). 

Fig. 2.— as function of R, for the three models (TCDM, ACDM, and OCDM). 

Fig. 3.— {{Sp/pY)r as function of R, for the three models (TCDM, ACDM, and OCDM). 

Fig. 4.— {61)r as function of R, for the three models (TCDM, ACDM, and OCDM). 

Fig. 5.— {v^)r as function of R, for the three models (TCDM, ACDM, and OCDM). 

Fig. 6. — P{6h\v) r, for ko = 0.01 /iMpc"^ and i? = 40 h~^Mpc. The solid and dashed lines are 
the distributions given by Eq.(^), with Xc = 10 {Ac = 0.1ai/a2) and Xc = I {Ac = 0.1/0.2) 
respectively; the dotted lines are Gaussian distributions with the same variance [given by 
Eq.(|6D] for Xc = 10. 

Fig. 7. — Same as Fig. 6, for R = 100/i~^Mpc. 
Fig. 8. — Same as Fig. 6, for R = 500/i~^Mpc. 
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Fig. 1. — Window functions of (solid line), {{Sp/pf')^ (dotted line), {5'q)r (short- 

dashed hne), and (v^)ij (long-dashed hne). 
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Fig. 2.— {5jj)R as function of R, for the three models (TCDM, ACDM, and OCDM). 
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Fig. 3.— {{5plpf)R as function of R, for the three models (TCDM, ACDM, and OCDM). 
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Fig. 4.— (51) R as function of R, for the three models (TCDM, ACDM, and OCDM). 
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Fig. 5.— (v^)^ as function of R, for the three models (TCDM, ACDM, and OCDM). 



- 26 - 



8 



4 - 



T I r 



T I r 



T I r 



TCDM n = 0.7 



ko=.01h/Mpc 



R=40Mpc/h 





-0.4 



J L 




^ I L 



-0.2 







0.2 



0.4 



H 



Fig. 6. — P{Sh\v) r for = 0.01 hMpc~^ and i? = 40 h~^Mpc. The solid and dashed hnes are 
the distributions given by Eq.(^), with Xc = 10 {Ac = 0.1ai/a2) and Xc = I {Ac = ai/a2) 
respectively; the dotted lines are Gaussian distributions with the same variance [given by 
Eq.(|6D] for Xc = 10. 
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Fig. 8.— Same as Fig.6, for R = 500 /i-^Mpc. 



